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^ ■ Abstract 



■ We investigate the relevance to the physics of under doped YBa2Cu306+x and 

■ YBa2Cu408 of the quantum critical point which occurs in a model of two an- 
tiferromagnetically coupled planes of antiferromagnetically correlated spins. 

m . 

0^ . We use a Schwinger boson mean field theory and a scaling analysis to obtain 



the phase diagram of the model and the temperature and frequency depen- 



. dence of various susceptibilities and relaxation rates. We distinguish between 

O ■ a low cj, T coupled-planes regime in which the optic spin excitations are frozen 



out and a high oj, T decoupled-planes regime in which the two planes fluctuate 
independently. In the coupled-planes regime the yttrium nuclear relaxation 
rate at low temperatures is larger relative to the copper and oxygen rates 
than would be naively expected in a model of uncor related planes. Avail- 
able data suggest that in YBa2Cu408 the crossover from the coupled to the 
decoupled planes regime occurs at T700A' or T 200K. The predicted corre- 
lation length is of order 6 lattice constants at T = 200^^. Experimental data 
related to the antiferromagnetic susceptibility of YBa2Cu408 may be made 
consistent with the theory, but available data for the uniform susceptibility 



1 



are inconsistent with the theory. 
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I. INTRODUCTION 



In this paper we investigate the relevance to the physics of underdoped YBa2Cu306+a; 
and YBa2Cu408 of a T = order- disorder transition which occurs in a model of two anti- 
ferromagnetically coupled planes of antiferromagnetically correlated spins. This transition 
might be relevant because: (a) in these compounds the basic structural unit is a pair of 
Cu02 planes separated from each other by the relatively inert CuO chains [ 0], (b) there is 
evidence for strong antiferromagnetic correlations within a Cu02 plane [0-^, (c) neutron 
scattering experiments find that for all uj and T studied a spin in one plane is perfectly 
anticorrelated with the nearest neighbor spin on the nearest-neighbor plane [ H,^, (d) for 
T < 150 K both the static uniform susceptibility and the various NMR relaxation rates drop 
rapidly as T decreases [ || suggesting [ 0-^] that the system is evolving as T is decreased 
towards a quantum disordered ground state with a gap to spin excitations and (e) the spin 
physics of La2_xSrxCu04, in which the Cu02 planes are very weakly coupled, is apparently 
rather different [ 0, suggesting that the behavior of YBa2Cu306+x may be due at least 
in part to a coupling between planes. A preliminary version of this work was published 
previously [ ] . 

The model we consider is a Heisenberg model of spins sitting on sites of the lattice 
depicted in fig. and has two coupling constants, both taken to be antiferromagnetic. One, 
Ji, couples nearest neighbor spins in the same plane. The other, J2, couples a spin in one 
plane to the nearest spin in the other plane. The Hamiltonian is 

h = j,y: st^-sf^j.Y.s^'^-sf'^ (1-1) 

<i,j>,a i 

Here i labels sites in a given plane, i and j are nearest neighbors in the same plane, and 
a = 1,2 labels two planes. There are two dimensionless parameters: J2/J1 and S, the 
magnitude of the spin. At temperature T = eq. ( |1 . 1|) has two phases. One is antifer- 
romagnetically ordered; the other is a singlet phase with a gap to excitations and no long 
range order. Varying J2/J1 and S produces transitions between the phases. In a physical 
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Heisenberg model S would take only half-integer or integer values 1/2, 1, 3/2, .... However, 
it is interesting to consider values S < 1/2 because the singlet phases occurring for small 
5* in this model may be useful representations of the "spin gap" behavior occurring in the 
underdoped YBa2Cu306+x- We shall determine the phase diagram and discuss the physical 
properties at the transition and on the disordered side. The T = transition is in the 
universality class of the 2+1 dimensional Heisenberg model and some universal properties 
have been determined [ |lO[|ll|. Here, we pay particular attention to the effects of the inter- 



plane coupling, J2. We show that one must distinguish between a low u, T "coupled planes" 
regime and a high uj,T "decoupled planes regime". In the coupled-planes regime, one linear 
combination (essentially the optic mode of the spin excitation spectrum) is frozen out and 
the low energy physics is determined by acoustic spin fluctuations in which moments in the 
two planes fluctuate coherently. In the decoupled planes regime the two planes fluctuate 
essentially independently. One important feature of the quantum critical point considered 
here is the close relationship between the susceptibility at small q and at a g near the or- 
dering wavevector [ |TT| . We show that this, when combined with interplanar coupling, has 
a surprising implication for the yttrium relaxation rate: in the coupled-planes regime it is 
larger, relative to the other rates, than one would expect from a model of uncoupled planes. 

We emphasize that eq. (|1 . 1|) is not a completely realistic model of YBa2Cu306.6 or 
YBa2Cu408 because it omits the itinerant carriers which make these materials metallic 
and indeed superconducting. The itinerant carriers also strongly affect the magnetism. In 
all of the hole-doped CUO2 compounds, long-range magnetic order disappears essentially at 
the metal- insulator transition [ |1^. It therefore seems likely that the itinerant carriers sub- 
stantially weaken the in-plane magnetic correlations. Itinerant carriers also presumably give 
rise to a particle-hole continuum of incoherent spin excitations which may strongly affect 
the physics and which in appropriate circumstances may change the universality class of the 
transition [ 0,|1^]. It has however recently been argued that there is some evidence that the 



behavior of YBa2Cu306.6 and YBa2Cu408 is in the universality class of eq. ( |1 . 1| ) [ 

The rest of the paper is organized as follows. In section II we solve eq. (|1.1| ) via the 



Schwinger-boson mean field theory. In section III we give the results of the Schwinger 
boson mean field theory for the phase diagram and the physical susceptibilities relevant to 
NMR and neutron scattering experiments. In section IV we combine selected results of the 
mean field theory with other arguments to map the problem onto a recently constructed 
scaling theory of the transition [ |TT| and give the relevant results of the scaling analysis. 



Section V is a conclusion in which we discuss the results and their relation to experiments 
on YBa2Cu306.6 and YBa2Cu408. It may be read independently of the previous sections by 
readers uninterested in the derivations of the results. Appendices present details of various 
calculations. 

II. MEAN FIELD SOLUTION 

We first study eq. ( p. . 1| ) by the Schwinger-boson mean- field method [ In this method 
one introduces bose operators which create a state of spin a on site i of plane a. One re- 
stricts oneself to the subspace in which each site on each plane has 25* bosons, corresponding 
to spin S] thus we enforce the constraint 

T.bli:^tJ = l + 2S (2.1) 

a 

A spin operator is written 

st^ = EblL^^^^^b,^,''^ (2-2) 



We now substitute eq. (|2.2| ) into eq. ( p..l|) . Then on the even sublattice of plane 1 and the 
odd sublattice of plane 2 we make the time reversal transformation which in the S* = 1/2 



case is: 



b\ - -b\ 

b\^ b\ (2.3) 



Rearranging and using eq. ( P?T| ) leads to 
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<i,jaal3 \ / \ 



(2.4) 



To proceed with the approximate mean-field treatment one introduces a Lagrange mul- 
tipher /i to enforce the constraint, and an in-plane bond field Q and a between- planes 
bond-field A to decouple the quartic interactions in eq. ( |1 . 1|) . In the mean field approxima- 
tion fi, A and Q are taken to be constant in space and time. The resulting theory may be 
diagonalized. The manipulations are standard and are given in Appendix A. The result is 
a model of two species of bosons (s, for symmetric under interchange of planes and a, for 
antisymmetric under interchange of planes) governed by the Lagrangian 

C'B = y] sl^[dr + UJk]Ska + [^^r + UJk+p]aka (2.5) 



ka 



with 



and 



^k = V/^' - (Q7fe + A)2, (2.6) 



7fc = ^(cos + cos ky) (2.7) 



(tt.vt). (2.8) 



Note that jk+p = —Ik- The parameters /i, Q and A are determined by the mean- field 
equations 

/^.coth^ = l + 2^, (2.9a) 



(27r)2cufc 2T 

£k W7. + A)7.^^^j^^^^/2J„ (2.9b) 



(27r)2 ujk 2T 
(Pk Q-fk + A ^, cufc 



(2.)^ coth- = 2A/J,. (2.9c) 

These equations are derived in Appendix A and imply Q,A,fi are real, and A, /i > 0. 
Further, one may change the sign of Q by shifting the origin of reciprocal space to k = P 
and interchanging the labels s and a. Thus we take Q > with no loss of generality. 
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The integrals in eqs. (|2.9|) depend on k only via 7^; one may therefore recast them as 
/ — > / d'yNij) where A^(7) is a density of states which is constant near the band edges 
7 = ±1 and logarithmically divergent at the band center 7 = 0. To obtain an analytically 
tractable model we replace A^(7) by 1/2. The resulting equations are solved in Appendix B. 
At T = we find the phase diagram shown in fig. 2. At J2 = we have two decoupled planes. 
As is well known [ |lO|,|lT],|l5|,|l6| a single plane has a transition dX S = S^ = {n /2- 1)/2 ^ 0.28 
(the numerical expression for Sc is obtained from the mean-field calculation) between a 
large-S* ordered state and a small-S" singlet state with a gap to all spin excitations. If one 
increases J2 in the ordered {S > Sc) phase of the one-plane model, one reaches at a J2 ~ Ji 
another transition line, at which the ordered phase is destroyed in favor of singlets which 
are principally between planes. At 5 = 1/2 the Schwinger boson mean field method yields 
a second order transition at J2 = 4.48 Ji. A previous series expansion study of eq. (|1.1| ) by 
Hida [ |17| yielded a second-order transition at J2 = 2.56 Ji and a very recent Monte-Carlo 
study by Sandvik [ |I^ found a second order transition at J2 = (2.7 ± 0.2) Ji. 

Interestingly, for S < Sc, we find that increasing J2 from J2 = initially reduces the 
gap in the singlet phase, thus moving the system closer to order. For S* < S < Sc (with 
S* = 0.19 in the mean field calculation) the phase diagram is reentrant. We believe that the 
physics behind the reentrance is that a small J2 splits the spectrum into acoustic and optic 
sectors, and because the optic sector involves a coupled motion of the spins in the two planes 
the effective spin of the model describing the low energy fluctuations is increased, promot- 
ing order, whereas at large J2 the between-planes interaction produced singlets, favoring 
destruction of order. 

In the disordered phase, /i > {Q + A) and the excitation spectrum has a gap at all 
wavevectors. Two gaps that will be particularly important in what follows are uj^ and 
given by 

00+ = _ (g + A)2 (2.10) 

and 
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= ^ix^ - (Q - A)2. (2.11) 

For the s-bosons uJj^ is the gap at A; = and is the gap at k = P; for the a-bosons and 
P are interchanged. We shall see that the matrix elements coupling the bosons to externally 
applied fields have a strong fc-dependence, so that measurable susceptibilities near k = 
differ dramatically from those near k = P. Both gap parameters are temperature dependent. 

At the T = phase transition, vanishes while uj_ > 0. Near the large J2 boundary 
of the ordered phase, we have A ~ Q and C(j_ ~ Ji ^ u^. The s-boson mode has one low 
energy branch, centered at k = 0, 

uj+{ky = iul + 2/i2(l - 7;.) =ujI + v^P (2.12) 

Here we have expanded 7^ = 1 — k'^/4:, set Q + A = /i and defined = /i^/2. 

In the lower, reentrant branch of the phase diagram, i.e. at T = 0, S* < S < Sc, J2 ^ Ji 
we find from eqs. (p7| ) that the phase boundary is given by 

^ = 7t{S, -S) + (tt^ - 8){S - S,)^ (2.13) 

On the phase boundary, 

= J2 + ... (2.14) 

Here the ellipsis denotes terms of order {Sc — 5*)^, {J2/ JiY and higher. If we tune through 
the phase transition by varying J2 we find that sufficiently deep in the ordered phase uj- 
increases as (J2-^i)^^^ as expected from spin-wave theory [|^, while uj_ approaches ojj^ very 
rapidly as J2 is decreased into the disordered phase. Indeed, within mean field theory we 
find find that for J2I Ji < 7r(S'c — S) — (8 — n'^/2){Sc — S)"^ a solution with cu^ 7^ c<j_ is 
not possible. The sharp transition from a solution with u!-uj+ to one with a;_ = uij^ is an 
artifact of mean field theory, but the qualitative result, that [u^ — uj+)/uj^ drops rapidly as 
one moves into the disordered phase, is likely to be correct for this model. The argument 
is that the ground state of the one plane model in the disordered phase is a singlet with a 
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gap of order Ji{Sc — S); for J2 less than this value the interplane coupling can only slightly 
perturb the singlets. 

For S = Sc and J2 ^ Ji the s-boson has two low energy branches, one centered at A; = 
with dispersion given by eq. (|2.12| ) and one given by 

u^kf = cu! + 2/i2(l + 7fc) = cu! + v^{k - Py (2.15) 

Here v = fi/2 + ... and the ellipsis indicates terms of order {J2/ Jif- 

We now consider T > 0. In a realistic model there are no phase transitions, however 
different regimes of behavior exist. In the mean field theory, crossovers between different 
regimes sometimes appear as unphysical phase transitions. We are interested in properties 
in the disordered regime near the critical line. At the T = phase transition, uoj^ vanishes; 
close to it is much smaller than Ji. For T < co'+(T = 0) the number of thermal excitations 
is negligible; the physics is of a singlet ground-state with a g-dependent gap to excitations. 
In the literature this is referred to as a "quantum disordered" regime. In this regime the 
low energy spectrum only involves one linear combination of the spin excitations in the 
two different planes; the antisymmetric one for k near P and the symmetric one for k near 
0. For ujj^{T = 0) < T < uj^(T = 0) this one linear combination becomes thermally 
excited. We refer to this as the "coupled-planes quantum critical regime". It only exists if 
uj^(T = 0) — cj+(T = 0) is large enough. In the coupled-planes critical regime uj^(T) ~ T, 
but uj^(T) takes its zero temperature value. The crossover from the quantum disordered 
regime to this quantum critical regime is identical to that occurring in a one-plane model. 
Finally, as T is increased through uj_ (T = 0) the other linear combination of spin excitations 
also becomes excited and the two planes begin to fluctuate more or less independently. If 
uj_{T = 0) ^ Ji then we find in the mean field theory that the behavior in this regime 
will be controlled by the T = critical point of a single plane. We refer to this regime 
as the " decoupled-planes critical regime". In the mean field theory the change from the 
coupled-planes to the decoupled-planes regime occurs via a second order phase transition. 
We expect fluctuations not included in the mean field theory will convert this into a smooth 
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crossover. In the decoupled-planes regime, both uj+ and uj- are proportional to T with the 
same coefficient. Fig. 3 shows the T-dependence of ujj^ and a;_ for parameters chosen so 
that uj+{T = 0) = 0, and indicates the different regimes. 



III. PHYSICAL QUANTITIES 

To obtain the magnetic susceptibihties we compute the hnear response of the system to 
an externally applied magnetic field /ij"'*. The details are given in Appendix C. We find it 
convenient to decompose the externally applied field into parts symmetric and antisymmetric 
under interchange of planes, and compute the linear response to 

M1)_M2) + 

Ag = E ' 2 ' 2 (3.1) 

Here and 0° are operators creating spin fluctuations symmetric and antisymmetric under 
interchange of planes respectively. The only non-zero susceptibilities are 

xr(^) = I rfte<'^+-)* < [Ol{t),0%m > (3.2a) 

and 

fOO 

Xl'ico) = cite^(-+^^)* < [Of (t), 0^,(0)] > (3.2b) 

The uniform susceptibility (written as a susceptibility per spin and not as a susceptibility 
per unit cell) is 

X{T) = hmx-iuj = 0) = ^Esinh-^ (^) (3-3) 

k 

We note that in the approximation of section II and the Appendices, namely. 



27r 

eq. (|373|) becomes 



kdk , , Ludu ^, 
N{-f)d-f = (3.4) 



roo rJ'T'T' poo r^TT 

X{T) = 2T[ f -1^] (3.5) 

Juij. 2T smh ix) Juj- 2T smh ix) 



'u)+/2T sinh (x) Ju)-/2T sinh (x) 
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From eq. (|3.5|) we see that in the coupled-plane quantum disordered regime {uj+(T = 0) ^ 
T), x(T) ~ u^e~'^+^'^] this is a factor of two smaller than would be found for a single 
plane in the disordered regime with the same gap. In the coupled-plane critical regime 
uj+(T = 0) ^ T ^ Lj^iT = 0), x(T) ~ T and in the decoupled-plane critical regime 
T ^ c<j_(T = 0), we also have x(T) ~ T but with a coefficient larger by a factor of two. 
This behavior is easy to understand: in the coupled-plane regime one of the two spin degrees 
of freedom per unit cell is frozen out; in the decoupled-plane region both are free to fluctuate. 
Next we consider a;) for q = P. From the results of eq. ( |3.2a| , p.2b| ) we flnd, for a;0 



xf'iu, T) = - coth(cu/4T) [Qiu - 2uj_) + Qiu - 2cu+)] (3.6) 

UJ 

Xp\uj, T) = —\Q{lo - K + c^_))[l + h{{uj'' ~ujI + ujI)/2uj) + h{{to^ + u^l - uj'i)/2uj)] 

UJ 

+ Q{u--uj+- u)) [h{{ujl -uj\- uj^)l2uj) - h{{ujl_ -uj\^ ^^)/2cu)]] (3.7) 

These formulae are plotted in flg. 5 for the parameters used to construct flg. 3. The 
sharp onset at T = is an unphysical feature of the dynamics in the mean fleld theory, 
which is removed when fluctuations are included [ |rT],|TU[. A related peculiarity of the 
mean fleld theory is that x'(g, = 0) decays as 1/g for q near P. We believe that despite 
the obvious artiflcialities the mean fleld expressions give some reliable information about 



the spin excitations of the model, as in the case of the one-plane Heisenberg model [ |T5 
In particular, we believe that the location of the peaks and their relative magnitudes and 
temperature dependences give a reasonable representation of the location, relative magnitude 
and temperature dependence of the peaks in the appropriate susceptiblities of the model. 
We see that from eqs. (|3.6|) and (|3.7|) and flg. 5 that at low T the important energy 
scale for Xp"" is 2ujj^ while for Xp " it is ujj^ ^ uj_. Everywhere in the disordered phase the 
difference between these energies is less than J 2- It is also clear that until the temperature 
becomes comparable to the scale T* at which the crossover from the coupled-planes to the 
decoupled-planes critical regime occurs, Xp " ^ Xp"; and that Xp" is conflned mostly to 
high frequencies, of order 2 J2. For this reason we suspect that the neutrons scattering data, 

11 



which have not observed any Xp" component for cu < AO meV, do not set a stringent hmit 
on J2, although they suggest that it is greater than 20 meV. 

We now turn to the relaxation rates. Details of the computations may be found in 
Appendix C. We begin with that of Yttrium, 1/^TiT, which is found to be 

^ = 2nD'\un-Y: 1^(^)1 V'M (3-8) 

The form factor g is equal to 4 for q near and to {q — P)^ for q near P; a more precise 
formula is given in eq. ( |C15|) . 

Using eqs. ( |3^ ) and ( |C15| ) we find 

dxx'^ /-oo dxx'^ 



+ 



(3.9) 



L+{T)/2T sinh'^{x) J<^„(T)/2T sinh'^{x) 

In the mean field approximation to both the coupled planes and the decoupled planes critical 
regimes, the yttrium rate 1/^TiT ~ T^, so it vanishes faster than the static susceptibility, 
while in the low T limit 1/^TiT ~ e~^+/^ is proportional to the uniform susceptibility. 
Note that the ratio of 1/^TiT to is the same in the coupled-plane critical regime as it 
is in the decoupled-plane critical regime. There are two compensating effects at work: in 
the coupled-plane regime there is only one spin mode at small q (instead of the two that 
occur in a model of two uncoupled planes) but in this mode the 8 spins add coherently to 
the yttrium relaxation rate, whereas in the decoupled-plane regime there are twice as many 
spin excitations (so one would naively expect the rate to be four times as large) but the two 
planes add incoherently (reducing the rate by a factor of two). 1/^TiT is plotted in fig. 4 
for the parameters used to construct fig. 3. 

We now consider the oxygen rate. This has one contribution from the small q fluctuations 
and one from fluctuations with q near P. The latter is suppressed by a form factor because 
each oxygen sits symmetrically between to copper sites, so a perfectly antiferromagnetic 
fluctuation would cancel on the oxygen site [ 0]. The result, derived in eq. (|C17|) is 



o^ = ^lim^E l/(9)r(xf M + xfM) (3-10) 
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This is precisely the usual result [ . Again using the approximation of eq. ( p.4| ) , we get 



^TiT Juj.(T)/2T sinh (x) Jc^_(t)/2t sinh (x) 

(3.11) 

In contrast to the yttrium, where the variation was by a factor of two, here the coefficient 
of the term varies by a factor of four between the coupled and decoupled-plane critical 
regimes. 

In the mean field theory the oxygen rate varies as in the critical regime. We shall see in 
the next section that in the more physically reasonable scaling analysis the antiferromagnetic 
spin fluctuations lead to an oxygen relaxation rate proportional to T in the critical regime. 

We now turn to the copper rate, which is given by 

c^yr = \^j^ol - 4i?7(g))'(x^^"(g, ^) + x'-^'^q, ^)) (3.12) 

Again this may be evaluated, giving 
1 {A-ABf " 



CuT^T 4 



/ dx/ sinh^(x) + 3 / dx/ sinh^(a; 

Juj4-(T)/2T Ju)-(T)/2T 



(3.13) 



/uj+(T)/2T Ju)-(T)/2T 

This also is plotted in fig. 4 for the parameters used in fig. 3, and again there is a factor of 
four change from the coupled-planes critical regime to the decoupled-planes critical regime. 
Note that in the critical regime the Cu relaxation rate is T-independent. Again, this is due 
to the artificiality of the Schwinger-boson mean field theory. The scaling analysis predicts a 
1/T behavior [ |ll[ 

IV. SCALING ANALYSIS 

The mean field analysis of the previous sections is known [ |ri|,|l^,^ to give incorrect 
results for dynamical susceptibilities for the single-plane Heisenberg model at finite spin 
degeneracy A^. A scaling theory which corrects these discrepancies has been developed [ 
TU|,|rT|. In this section we extend the scaling theory to the coupled-plane system of interest 



here, and also construct a universal amplitude ratio for the NMR Ti and T2 relaxation times. 
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In the scaling theory of the one-plane model the important parameters for the physics 
on the disordered side of the phase boundary are the spin- wave velocity, v, the T = gap to 
spin-1 excitations, A, the quasiparticle residue of the lowest-lying S=l excitation at T = 0, 
Aqp and the temperature, T. Low energy physical quantities are universal functions of these 
parameters. Further, one must distinguish the "quantum disordered" T <^ A regime from 
the "quantum critical" T ^ A regime. 

We now extend the theory to two coupled planes. We expect on general grounds, and 
showed explicitly using the mean field theory, that the between-planes coupling J2 splits 
the spin excitation spectrum into acoustic and optic branches. The minimum gap to optic 
excitations, Aq, is nonzero for J2 > 0; the acoustic excitations acquire a gap A^ in the 
disordered phase. In the mean field treatment of the previous section, A^ = 2co'_|_ while 
Ao = u!+ + uj_. We found that on the phase boundary Ao ~ J2 while A^ = 0. As one moves 
into the disordered phase, A^ rapidly approaches Aq. Deep in the ordered phase, for small 
J2, Ao ~ y/JoJi in agreement with spin- wave theory [ |[ 

As we tune the system through the order-disorder transition, Aq remains nonzero, so the 
optic fluctuations are irrelevant in the renormalization group sense. Within the mean field 
approximation the transition is thus one in which the Heisenberg order vanishes and a single 
two-fold degenerate spin- wave mode of velocity Va acquires a gap A^. The usual universality 
arguments then imply that the transition is in the previously considered universality class, 
and that at u;, T ^ Ao the physical quantities are given by the previously calculated univer- 
sal functions [ |ll| evaluated at arguments Va and T/A^. We refer to this ti;,T ^ Aq regime 
as the coupled-plane regime, and distinguish the T ^ Aq ^^coupled-plane disordered regime" 
from the T ^ A^ ^^coupled-plane critical regime" . Of course, if (Ao — Aa)/Ao is too small 
the coupled-planes critical regime may not exist. 

There is one subtlety in the analysis of the coupled-planes regime: the universal forms 
give susceptibilities in units of emu/ area. A unit cell contains two Cu atoms (one in each 
plane); the susceptibilities are evenly divided between the two planes, thus the susceptibilities 
per Cu are one half of the one-plane values. Further, because the Cu and O nuclei of interest 
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for high-Tc NMR experiments sit in one Cu02 plane, the hyperfine coupling of the Cu or O 
nucleus to the surviving spin-degree of freedom is half as large as in a single-plane model, so 
the Cu and O relaxation rates, which go as the square of the hyperfine coupling constant, will 
be one quarter of the usual size. The yttrium nucleus, however, sits between two planes and 
over the center of a plaquette of four Cu nuclei. In the coupled-planes regime the hyperfine 
coupling per moment is half the expected value, but the eight nearest neighbor spins add 
coherently , so the yttrium NMR rate is only one half of the expected value. These factors 
were explicitly derived in the mean field analysis of the previous section; we have given a 
qualitative argument here. 

As one increases the temperature from zero in the disordered phase one passes first 
through the coupled planes disordered regime and then, if (Aq — A„)/Ao is sufficiently large, 
through the coupled-planes critical regime. As one continues to increase the temperature 
it becomes comparable to the optic mode energy Aq and the bilayers become uncorrelated. 
If Ao is sufficiently small, i.e. if J2 ^ Ji then at TA^ the physics will still be controlled 
by a T = critical point. We showed using the mean field theory that in this regime 
each plane fluctuates independently and is in the quantum critical regime of a one-plane 
model. We believe this conclusion survives beyond mean fleld theory. We refer to the 
regime TAq as the decoupled-planes critical regime" . It is probable that one could extend 
the calculation of reference [ |lT| to incorporate fluctuations into our mean- fleld analysis of the 
crossover between the coupled and decoupled planes regimes, but we have not attempted 
this. Instead, we consider the uniform susceptibility and NMR in each plane separately. 
The scaling properties of the single-plane transition have recently been elegantly derived 



and discussed [|TT|. We summarize, as briefly as possible, the results we will need and their 
extension to the two-plane system. 

The correlation length , ^ is a universal function of T and A, 

e = ^X{ksT/A) (4.1) 

with X{y) a universal function which tends a number very nearly unity (1.03 in a large N 
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expansion with terms of order 1 and 1/N included and N set equal to 3) as y tends to infinity 
(so ^ ~ hv/ksT for T ^ A) and to y as y tends to zero (so ^ tends to hv/A for small T). 
Note that for J2 ^ Ji, v is essentially the same for acoustic and optic modes, so that this 
relation is valid in both coupled and decoupled-plane regimes. 

The divergent part of the order parameter susceptibility, Xaf, may be written 

XAFiq, co; T) = XAFe-yAPiq^, cu^/v; T/A) (4.2) 

Here xaf contains the dimensions and (pAF is a universal function. We have introduced the 
exponent r] for completeness even though for the transition in question it is very nearly zero 



[ |ri|,^. In the coupled-planes regime we should interpret this as a susceptibility per unit 
area; in the decoupled-planes regime as a susceptibility per area per plane. Here xaf is a 
critical amplitude; its value is not universal but from it and other measurable quantities 
universal amplitude ratios may be constructed. In particular, xaf = Agp/v^ where Agp is 
the quasiparticle residue of the lowest-lying S=l excitation in the disordered phase at T=0 
and we have made explicit the factors of the velocity which were set to unity in [ |rT| . 

The uniform susceptibility has also a scaling form; here even the dimensional prefactor 



is known [ |TT|. The susceptibility in units of emu per area is 



Here g is the electron g-factor and /is is the Bohr magneton. Again this is a susceptibility 
per unit area in the coupled planes regime and a susceptibility per area per plane in the 
decoupled planes regime. Dg is the spin diffusion coefficient; it has a dependence on T and 
q which we have suppressed here. 0„ is a universal function which tends to a constant as 
y tends to infinity and vanishes rapidly as y tends to zero. It is usually assumed that the 
quantity g^B niay be taken to have the free electron values because in the nonlinear sigma 
model treatment of the critical point it is not renormalized from its bare value; however this 
assumption has not been proven for more general models. 

We shall be most interested in applying this formula in the critical regime T ^ A; we 
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therefore proceed to estimate Ds in this regime. We argue that for the modes relevant to 
NMR experiments, 

= Dy (4.4) 

Here -D^ is a number, presumably of order unity. The factor of 1/T arises as follows. 
For excitations of velocity v and scattering rate F, Dg ~ v'^/T. Further, conventional 
dynamic scaling suggests that modes at scales shorter than the correlation length are weakly 
damped, while those at longer scales are overdamped, implying that the scattering rate for 
spin excitations is proportional to T. The factor of In^^^ comes from the breakdown of 



hydrodynamics in two spatial dimensions [ ^ 



We are now able to discuss relaxation rates. We begin with the yttrium rate. We saw in 
the previous section that the yttrium nucleus is coupled only to fluctuations of the uniform 
magnetization. By comparing the notations of this and the previous section we see that in 
the coupled-planes regime the coupling constant is ADa'^/fisg (recall that there are 8 Cu 
neighbors but that the spin density is evenly divided between the two planes). Calculating 
the relaxation rate in the usual way gives 

y^, ^ = .3no zc In'/' IM 4.5 

^t^kbT h Dy$, IT 

Recall that refers to two planes; normalizing per plane restores the factor of 32 from 
3.9||. A power of In^ was obtained previously by Chakravarty and Orbach in a calculation 



of relaxation in an ordered magnet [ the different power comes because they did not 
consider the corrections to hydrodynamics. The logarithm will be cut off at some scale 
by a three dimensional coupling and is presumably not important in practice. The 
important result is that in the quantum disordered regime both Dg and go as 1/T so the 
yttrium relaxation rate in this model is proportional to up to logarithms, as was found 
in the mean field theory of the previous section. In the decoupled-planes critical regime the 
calculation is identical except that we must add two contributions, one from each plane. 
For each contribution the coupling constant is ADa? / and we must neglect interplanar 
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correlations. The result is a factor of two increase in the coefficient of the term, as was 
found in mean field theory. 

We now consider the oxygen relaxation rate. This has one contribution from the small-g 
fiuctuations which may be evaluated as we did for yttrium and which will be seen to be 
sub-dominant, and another contribution from the antiferromagnetic fiuctuations, which we 
evaluate. The coupling constant connecting the oxygen nucleus to an antiferromagnetic 
fiuctuation in a given plane of wavevector q (measured from P) may be written Cafcl^ fiQCi) 
where Caf is a priori not the same as the coupling constant C introduced before. Symmetry 
implies that f{x) ~ x"^ at small x. Combining this with eq. (|4.2| ) for xaf gives, in the coupled 
planes regime: 

h l^^^^'^^^.(e/a)-(^+'')0o(T/A.) (4.6) 



OTiksT 4 hv fil 
Here (j)o{z) is a universal function obtained by integrating limy^o x'^ (l)'^p{x, y; z)/y over x. 
We have assumed the integral converges; in the spin-only model this is reasonable because at 
momentum scales larger than the inverse correlation length the model goes over to spin-wave 
theory and the Integra Is there converge. In a more general model the issue of convergence 
is less clear. Thus the oxygen relaxation rate in this model scales at T^+'' with a non- 
universal prefactor involving both the hyperfine coupling Caf and the amplitude Xaf- The 
oxygen relaxation rate scales differently from the yttrium because the oxygen is coupled 
(albeit weakly) to the antiferromagnetic fiuctuations, while the yttrium is not. As we have 
previously argued, the constant Caf is larger by a factor of two in the decoupled-planes 
regime than it is in the coupled-planes regime, leading to a factor-of-four change in the 
relaxation rate. 

We finally consider the copper. This is coupled to the spins by a matrix element 
^afc^^/z^b- The relaxation rate in the coupled-planes regime is 

^-^"irii/n.im) (4.7) 

Thus the Cu relaxation rate 1/TiT in this model scales at T^~^ times a non-universal 



prefactor involving both the hyperfine coupling and the amplitude Xaf- The formula p7 



has been previously given [ 0. The same factor of four change in the coefficient of the 
leading T-dependent term between the coupled-plane and decoupled-plane regimes that 
occurred for the oxygen relaxation rate occurs for the copper. 

The Ti relaxation rate is determined by the imaginary part of x- The T2 rate measures 
the real part of x- Specifically, in circumstances relevant to experiments on high-Tc materials 



23 



d)^ = n„5: [AAFX{q.^ = Q)? (4.8) 

^2 q 

where rim is the density of NMR nuclei. Substituting the scaling ansatz and integrating 
gives 



rp ""m 2 



2 (e/«)^-VT2(T/A) (4.9) 



By combining eqs. ( [4.7| , [4.9| ) we see that apart from the factor and the velocity f , 
the ratio of T2 to TiT is universal, and indeed takes the same value in the coupled-planes 
and decoupled-planes critical regimes. Note however that in the coupled-planes regime the 
contribution of the optic excitations to I/T2 will be large (of order I/Aq). The contribution 
of the acoustic sector is of order I/Aq. Thus I/T2 will attain its universal value in the 
coupled-planes regime only if (T, A^) <^ Ao 

Sokol and Pines [ ^ have previously made the interesting observation that the observed 
T-independence of T2/T1T in YBa2Cu306.6 for T > 150 K suggests that the magnetic 
dynamics in this material is controlled by the z = 1 critical point considered here. We 
see that the magnitude provides information about the velocity, v, and that this must be 
consistent with the uniform susceptibihty. 



V. CONCLUSION 

We have studied some aspects of the T = magnetic-non-magnetic transition occurring 
in a model of two antiferromagnetically coupled planes of antiferromagnetically correlated 
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spins. The model is defined in eq. ( |1 . 1| ) and depicted in fig. 1. The two dimensionless 
parameters are J2/J1 (the ratio of the between-planes couphng J2 to the in-plane coupling 
Ji) and S*, the magnitude of the spin in one plane. The important dimensional parameter 
is the spin wave velocity, v. The phase diagram at T = in the J2/J1 - S plane is shown 
in fig. 2. For a single plane (i.e. J2 = 0), decreasing S through a critical value Sc causes 
a phase transition between a magnetically ordered phase and a singlet phase with a gap 
to spin excitations. Some properties of this transition were determined by Chakravarty, 
Halperin and Nelson [ ^ and it was analyzed in detail by Sachdev, Chubukov and Ye[ p!l| . 
For the coupled-plane system we found using a Schwinger-boson mean field theory that a 
large value of the interplanar coupling J2 destroys the magnetism even for SSc (because it 
favors binding of nearest neighbor spins on different planes into singlets), while a small J2 
promotes order by increasing the effective size of the spin in a unit cell. The interplay of 
these two different sorts of physics leads to the reentrant phase diagram shown in fig. 2. 
This phase diagram differs in an important respect from our previous interpretation of the 
data on spin susceptibilites of La2_xSr^Cu04 and YBa2Cu306.6 [ both compounds 

it is clear that doping destroys the magnetism. In La2_xSr^Cu04 interplane coupling is 
negligible and at least the Cu relaxation rate and single crystal susceptibility data show no 
clear evidence of a singlet phase with a gap to excitations (although some susceptibility and 
Knight shift data have been so interpreted [ p|,pT|). On the other hand, in YBa2Cu307_5 it is 
clear that the nearest-neighbor CUO2 planes are coupled by an interaction at least of order 
300 K, and "spin-gap" effects are very easily observed in susceptibilities and relaxation rates 
for 0.1 < 5 < 0.5. Thus it appears that in the real materials a presumably modest between- 
planes coupling promotes "spin-gap" behavior, and therefore that the spin-only model is 
missing some essential feature of the physics, most likely related to the presence of mobile 
holes. 

Although it is not completely realistic, the coupled-plane model might capture some 
aspects of the physics of YBa2Cu306.6 and YBa2Cu408. We therefore calculated the pre- 
dictions of the model for the temperature and frequency dependence of the susceptibilities 
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measured in NMR and neutron scattering using the Schwinger boson mean field method 
and a scahng analysis. We studied parameters such that the model has no long range order 
at T = 0. The behavior in the disordered phase of a single plane of Heisenberg spins is 



understood [ |TT],[T6|. In the single-plane case, the important parameter is the T = gap 
to spin one excitations, A. The spin-one excitations are essentially spin waves with a gap. 
There are several regimes of temperature. For T <^ A the number of thermal spin exci- 
tations is negligible; the static spin susceptibility at g = and the dissipative part of the 
dynamic susceptibility at all q and at <^ A have an activated temperature dependence 
~ g-A/T 'jj-^jg regime is referred to as the "quantum disordered regime". If the microscopic 
exchange constant J ^ A then for A ^ T ^ J another regime exists in which the physics 
is dominated by the T = critical point but the gap is not important. In this regime 
the static uniform susceptibility is proportional to T and the antiferromagnetic correlation 
length grows as 1/T. The regime is referred to as the "quantum critical regime". 

In the two plane model of interest here the between-planes coupling J2 splits the spin 
excitation spectrum into acoustic and optic modes. There are two important scales: A^, 
the T = gap to acoustic excitations and Aq, the T = gap to optic excitations. Both 
gaps are nonzero in the disordered phase. At the antiferromagnetic-singlet transition A^ 
vanishes while Aq remains non-zero. At the transition we found from the mean field theory 
that Ao = ^2- As one moves deeper into the disordered phase, (Aq — Aa)/Ao decreases 
rapidly. If one is sufficiently close to the phase boundary, so that < A^ ^ A^ ^ Ji, 
there are three regimes. These are depicted in fig. 6. For (T, u) <^ Aq the optic mode of 
the two plane system is frozen out and the physics is dominated by the acoustic mode of 
the two plane system. We refer to this as the " coupled-planes" regime. For (T, u) ^ A^ 
even the acoustic mode is frozen out. This is the "coupled-planes disordered regime". For 
A(j ^ T ^ Ao the acoustic mode is thermally activated and the system is in the "coupled- 
planes quantum critical regime". Finally, for Aq <^ T the coupling between the planes 
become negligible and the planes fiuctuate independently. This is the "decoupled planes" 
regime. If Aq ^ T ^ Ji then we argued that the spin dynamics in each plane is separately 
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given by the quantum critical behavior of a one-plane model. Of course if (Aq — Aa)/Ao is 
too small, the coupled-plane critical regime does not exist, while if Aq/ Ji is too large, the 
physics in the decoupled-planes regime will not be controlled by a T = critical point. 

We now summarize the results we obtained in sections III and IV for the temperature and 
frequency dependence of the NMR rates and susceptibilities. We found that the difference, 
Afj — Aa < J2 everywhere in the disordered phase, and that the contribution from the 
optic modes near the antiferromagnetic point was rather small and only weakly temperature 
dependent if Ao ^ A^. The uniform susceptibility, x(T), is activated in the coupled planes 
disordered regime; in the coupled-planes critical regime x{T) = 0.5ET (here E is a. number) 
and in the decoupled-planes critical regime x{T) = ET. The factor of two change in the 
coefficient of T between the coupled and decoupled planes regimes was derived explicitly 
from the Schwinger boson mean field theory. The physical origin, we believe, is that in 
the coupled-plane regime one of the two spin modes at each k (namely, the optic mode) is 
frozen out and does not contribute to x, whereas in the decoupled-planes regime both modes 
contribute. 

The T-dependences of the nuclear relaxation rates are more subtle. They are summarized 
in fig 6. We consider first the copper and oxygen rates. In the quantum disordered regime 
all rates are activated. In the coupled-planes critical regime the Cu rate 1/TiT = 0.25A/T 
while the oxygen rate is given by 0.25CT Here we have set the exponent rj, which is in 



practice very small [ ^ to zero. A and C are constants. In the decoupled planes critical 
regime the formulae are the same except that the number 0.25 becomes 1. The factor of four 
change in the coefficient of the leading T-dependence of the Cu and O rates between the 
coupled-planes and decoupled-planes critical regimes was derived from the mean field theory. 
We believe that the physical origin is that in the coupled-planes regime only the acoustic 
mode contributes to relaxation rates. In this mode the spin fiuctuation is even divided 
between the two planes; the hyperfine coupling constant connecting a Cu or O nucleus in a 
given plane to the spin fiuctuation is thus half of what it would be in a one-plane theory, 
and the rate goes as the square of the hyperfine coupling. The factor of four variation has an 
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interesting implication for estimates of the relative strengths of the antiferromagnetic spin 
fluctuations in different materials. Authors (including us) who had considered the question 
previously argued that because the Cu relaxation rate in YBa2Cu307_5 was smaller than 
in La2-xSrxCu04, while the hyperfine couplings were approximately the same, the spin 
fluctuations must be weaker in the former material. We see now that until one knows 
whether YBa2 Cus Oj^s is in the coupled-planes or decoupled-planes regime one cannot 
meaningfully compare the magnitudes of the relaxation rates to those of La2Cu04^. 

The yttrium relaxation relaxation rate behaves slightly differently, because the yttrium 
nucleus is coupled to a symmetric combination of Cu nuclei in the two planes. In the 
disordered regime the yttrium rate is activated, in the coupled-planes critical regime it goes 
as 0.5DT^ and in the decoupled-planes critical regime it goes as DT^. The yttrium rate 
changes only by a factor of two between the coupled-planes and decoupled-planes critical 
regimes because in the coupled-planes regime the spins in different planes move coherently 
while in the decoupled-planes regime the spins move incoherently. The difference between 
coherent and incoherent addition of spin fluctuations produces a factor of two which partially 
compensates for the factor of four discussed previously. 

We now consider the implications of our results for experiments on YBa2Cu306.6 and 
YBa2Cu408. Below T = T* ~ 150 K all of the relaxation rates including the Cu 1/TiT drop 
as T is decreased. We believe that this can only occur if the physics below T ~ 150 K is 
dominated by thermal excitations above a T = singlet state with a gap to spin excitations. 
In a spin-only model such as that considered here one would model this by choosing a value 
of spin S such that the system was near to, but on the disordered side of, the phase boundary 
in fig. 2. Further, the value of T* implies that ~ 150 K. We must next consider the 
value of the interplanar coupling J2. The experimental evidence is not conclusive. Only 
the acoustic excitation of the two-plane system has been observed via neutron scattering in 
any member of the YBa2 family of high-Tc materials [ The main effort has been at 
energies less than 40 meV and temperatures less than 150 K. This would suggest that in the 
metallic YBa materials, Aq > 40 meV. NMR measurements on YBa2Cu408 provide more 
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information. The Cu Ti has been found to obey very well the Curie law 1/TiT ~ 1/T for 



200 K < T < 700 K [ |2^. This implies that the temperature scale at which the bilayers 
become decoupled in this material is less than 200 K or greater than 700 K. We believe a 
bilayer coupling greater than 700 K would be very hard to justify on theoretical grounds. 
The argument is that the in-plane J ~ 1500 K, while J2 ^ Ji because band structure [ 
25| and photoemission [ |26| results imply that the between-planes hybridization is much 
less than the in-plane hybridization, and the exchange energy scales as a high power of the 
hybridization. However, a bilayer coupling much less than 200 K may not be consistent 
with the neutron data. From the mean field theory we found Aq < J2 + A^; our estimate 
Aa ~ 150 K then implies < 350 K if J2 < 200 K. Thus presently available data provide 
somewhat contradictory answers to the question whether YBa2Cu408 is in the coupled- 
planes or decoupled-planes regime for 200 K < T < 700 K. In what follows we consider 
both possibilities. 

We now turn to a more quantitative discussion. Sokol and Pines have proposed that the 
magnetic dynamics of underdoped high-T^ materials are determined by the critical point 
discussed here, and their discussion has been amplified and extended by Barzykin, Pines, 
Sokol and Thelen [ |^ . They do not consider the bilayer coupling. They propose that these 
materials are in the quantum critical regime for T > 150 K and in the quantum disordered 
regime for T < 150 K. The essential piece of evidence they cite in favor of their proposal 
is the observed approximate T-independence of the ratio T1T/T2 for 150 K < T < 300 K. 



The observed [ 27] magnitude of this ratio is approximately one-half of the value observed 



[ |2^ in La2Cu04 and calculated [ |29| for the S = 1/2 Heisenberg model with J = 0.13 eV. 
The density Um of ^^Cu NMR ions is the same in the La-Sr and YBa materials, so we 
conclude from eq. ( |4.9| ) that if the magnetic dynamics of YBa2Cu306.6 and YBa2Cu408 are 
well described by the universal scaling forms in either the coupled-planes or the decoupled- 
planes regimes, then the appropriate spin-wave velocity v is about one half of the 0.8 eV A 
appropriate for La2Cu04 [ .i.e. 
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WYBa = 0.4 [eV A] (5.1) 
The value of v makes a prediction for the magnitude of the correlation length. From eq. 



4.1|| we find ,^/a = 4atT = 300i^' and ,^/a = 6atT = 200/^; below this temperature the 
crossover to the quantum disordered regime presumably means that the T dependence of 
the correlation length becomes much weaker. These lengths are rather larger than observed 
in neutron scattering experiments [ . 

These considerations also have an implication for the Cu Ti relaxation rate. In the 
quantum critical regime, the theoretical result for the Cu relaxation rate is 1/TiT = A/T. 
The coefficient A ~ Agp/ [ 0] where Agp is the quasiparticle residue of the lowest lying 



S=l state above the gap at T = 0. It may in principle be determined from neutron scattering 
mea surements at low T. The coefficient A also depends on hyperfine couplings (which have 
been claimed to be the same in La2Cu04 as in YBa2Cu306.6 and YBa2Cu408 [ PI,|3T|) and on 
whether the material is in the coupled-planes or decoupled-planes regimes. In the absence 
of a measurement of Agp one cannot definitively calculate the Cu relaxation rate; however, 
it is interesting to attempt to estimate it. We first argue as follows: La2CuOi has been 
clai med to be in the quantum critical regime for T > 600K [ further, at these 

temperatures the Cu Ti has only a weak doping dependence, implying that the combination 
Agp/v"^ depends only weakly on doping. If this is c orrect, then we would expect that 
if YBa2Cu306.6 and YBa2Cu408 were in the decoupled planes regime, A would be well 
approximated by the value A = 3300 sec~^ [ ^ appropriate to La2Cu04, while if it were 
in the coupled plane regime we would expect A = 800 sec~^. (In fact, the values might be 
slightly larger since the J for YBa2Cu30Q,o is smaller than the J for La2Cu04 [ ^). In 
fact, in YBa2Cu408, A = 1600 sec~^ [ ^ compatible with neit her estimate. An alternative 
argument would be to say that Agp has the dimension of energy, and the characteristic 
energy scale is given by dividing the velocity, v, by the lattice constant, a, implying A ~ 1/u 



so from eq [5.1| we w ould expect that in the coupled- planes regime the value of the Cu Ti 



would come out about right. A third possibility is that for some unknown reason Agp could 
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drop by a factor of 8 upon going from La2-xSrxCu04 to YBa2Cu408, so tha t t he data 
for T > ISOi^ would be consistent with the decoupled planes regime. Of course a fourth 
possiblity is that the theory is not applicable. Neutron measurements of absolute scattering 
intensities on underdoped YBa materials would be very helpful in resolving this issue. 

Our discussion so far has emphasized properties related to the antiferromagnetic fluctu- 
ations. These are in a sense robust and relatively model- independent, depending as they 
do primarily on the existence of a growing correlation length, a weak J2 and a spin-gap at 
low T. We now consider small-g properties. In the Heisenberg model, the behavior of the 
small-g susceptibility is very closely tied to the behavior of the large-g susceptibility [ p!T| , |32| . 
In a more fermi-liquid-like model this need not be true, and therefore the bilayer coupling, 
which strongly affects the behavior near the antiferromagnetic point, need not in a more 
realistic model also strongly affect the susceptibility near q = 0. It is also clear that the 
spin-only model is not, by itself, a reasonable description of the small-g spin dynamics of 
Y Ba2Cu^0'j^s or YBa2Cu40s- For example, the observed yttrium nuclear relaxation rate 
does not vary as but is more nearly proportional to the static susceptibility [ Also, 
our calculated oxygen nuclear relaxation rate is too small by about a factor of 16 to explain 
the data [13^. Thus, we must invoke an extra contribution to x" existing at least at small 



q. It is natural to suppose that this is due to the mobile carriers and that thereofre the 
contribution to x" is of more or less the fermi-liquid form. However, the fermi-liquid like 
contribution cannot be appreciable near the antiferromagnetic point or it would overdamp 
the spin waves and change the universality class of the magnetic fluctuations [ [l^JT^. It 
has very recently been argued [ ^ that in the Shraiman-Siggia model of doped antiferro- 
magnets [ ?] precisely the required behavior occurs, with the fermions making an additive 
contribution to the small q but not the large q susceptiblities. A similar conclusion was 
drawn from high temperature series expansions [ ^ . In the resulting tw o-component pic- 
ture X = Xspin + Xqpj with the fermi-liquid like piece Xgp providing the yttrium and oxygen 
relaxation and the Xspin given by the theory we have discussed and providing the nontrivial 
temperature dependence of the uniform susceptibility and the Cu relaxation rate. Sokol 
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and co-workers have made a similar argument on phenomenological grounds, claiming that 
the temperature dependence of x YBa2Cu408 for 200 K < T < 700 K is consistent 
with the behavior of the spin-only model in the quantum critical regime [ Q . Now whether 
one supposes the material to be in the coupled-planes or decoupled-planes critical regime, 
the theory implies x{T) = a + hT . We ignore the value of a, on the grounds that it is 
dominated by the fermions which are beyond the scope of the theory, and consider the 
value of 6, which is given by eq. (|4.3|) . From eq. (??) assuming g^s takes the same value 
as in La2CuOi we find 6 = 5 x 10~^ states/eV — Cu — K for the decoupled-planes regime 
and b = 2.5 x 10"'^ states/eV — Cu — K for the coupled-planes regime. The data say that 
6 ~ 1.6 X 10^^ states/eV — Cu — K [ This too weak T dependence of x in the critical 



regime suggests to us that the straightforward two-component approach is not applicable to 
YBa2Cu408, and that the presence of carriers modifies the magnetic behavior more dramat- 
ically. We note, however, that two results which seem to be qualitatively consistent with 
the calculations presented here are: (a) although the difference is not dramatic, the yttrium 
rate seems to drop faster than the oxygen rate as T is lowered and (b) the magnitude of the 
yttrium rate is larger than expected from a model in which the planes are uncoupled [ . 

Whatever is the correct theory, the distinction we have drawn between the coupled-planes 
and the decoupled-planes regimes will still be important. Further, if the NMR and neutron 
data on underdoped YBa superconductors are both taken at face value, then t he crossover 
between the coupled-planes and decoupled-planes regimes occurs either at T > 700-?^ or at 
T ~ 2007^. The larger value seems to us to require an implausibly large between-planes 
coupling; the smaller value would imply that the crossover to the quantum disordered regime 
is complicated by a simultaneous freezing out of the optic mode of the bilayer system. 
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APPENDIX A: DERIVATION OF MEAN FIELD EQUATIONS 



We begin from eq. ( |1.1| ). We write the model as a functional integral, introduce a field 
Q<i,j> decouple the Ji interaction, a field Aj to decouple the J2 interaction and a field fj, 
to enforce the constraint. We find for the partition function, Z, 



Z = J VA+AVQ+QVb^bVi^exp J^^ drC'^ 



(Al) 



with 



+ 1 E e^&J?QS.> + h.c. 

<i,j>,aa 

+ EeeA. +h.c. 

ia 

+ E ^%^+^^ (A2) 

<i,j>a i ■J'i 

The different normalizations of Q and A have been introduced, so that the final expression 
for the boson energy, eq. ( |A9| ) , has no numerical factors. We next introduce symmetric (s) 
and antisymmetric (a) bose fields via 

^ia = ~'~ '^''^"^ (A3a) 

C^ = ^(4a-4a) (A3b) 

We make the mean field approximation of space and time independent Q,A,fi, Fourier 
transform the boson operators and obtain 

= E 4a[dr + /i]4. + l^iQlk + A)sls'l,, + h.c] 
ka 
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+ E 4Adr + /i]aL + [\{Qlk - A)aLa'i,„ + h.c] (A4) 

ka 

with 

The boson part of this equation may be decoupled by a Bogohubov transformation. The 
resulting quasiparticles s and a are defined by 

4a = cosh 6*^40 - ^Vah9kS-k,a 

4a = cosh 9k+pa\^ + sinh Ok+pa^k^a ( A6) 

with 

tanh2^fc = [g7fc + A]//x (A7) 

and 

P=(7r,7r) (A8) 

The energy of the s-bosons is 



ook = V/"' - {Qlk + A)2 (A9) 

The energy of an a-boson at a wavevector k is uJk+p- 

The free energy F may be computed in the standard way and is 

F = 4ArT^ln[2sinh(cjfe/2T)] + NQ"^ /2Ji + 2iVAV^2 - 2A^(1 + 25)/i (AlO) 

k 



The mean field equations, eqs. (|2.9| ), follow from differentiating this equation with respect 
to fi, Q and A. 

APPENDIX B: APPROXIMATE SOLUTION OF MEAN FIELD EQUATIONS 

We begin with eqs. ( p^) . We recast them as / — > / d'jN{'y), we replace 
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A^(7) by 1/2, we normalize Q and A by /i and integrate, obtaining at T = 

sin-i(A + Q) -sin"i(A-Q) 



sin-^(A + Q) - sin-^(A - g) + (A - Q)^! - (A + Qf - (A + g)^l - (A - QY _ 

4g2 ~ 2Ji 



2Q l+2S(Bla) 



(Bib) 
(Blc) 



l_(A_g)2_^l_(A + g)2 

2Q ~ J2 

These three equations may be reduced to one by taking the sine of eq. ( pia]) and substituting 



into eq. ( |Blb|) to obtain an equation for fi{Q), solving eq. (|Blc| ) to obtain an equation for 



A(g) and then substituting the results into eq. ( |Blb|) . However, for our purposes a simpler 
approach suffices. We first locate the critical point at which the minimum boson energy 
vanishes. In the notation of this appendix this implies A* + g* = 1 (we denote by * the 
values of the quantities at the critical point). Then eq. ( |Bla| ) may be solved for Q*. For 
S > Sc = {tt/2 — l)/2 this has only one solution; e.g. at S' = 1/2 

Q* ^ 0.277 (B2) 

For S* < S < Sc, with S* = 0.19 there are two solutions, one at Q near 1 which is the lower 
energy solution for Ji ^ J2 and one at Q near 1/2 which is the lower energy solution for Ji 
near J2/4. For S < S* there are no solutions. 

Once a solution for Q* is found, eq. ( [Blc| ) implies 

- ' (B3) 



J2 2^g*(i-g*) 

and eq. (pib|) implies 



J* Jl + 2S)^Q*il-Q*)-l + Q* 

Tt = ^ g^^ ^^^^ 

Solving eq. (|Bla| ) and then using the solution in eq. ( P^ ) yields the phase diagram given in 
fig. 2. 

We now consider the TO behavior. We are most interested in the regime near the phase 
boundary, and in small J2. We therefore solve the equations perturbatively in the small 
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parameters Sc — S, J2/ Ji and T. We neglect terms of third order and higher in these small 
parameters. To this order the T dependent terms may be evaluated exactly. The equations 
are conveniently expressed in terms of the variables uj+, uj^ and yU and are (note we need ix 
only to first order in the small parameters) 

TT - cu+//i - uJ_/^^ + {2T/^^){f{uJ^/T) + f{uj_/T)) , 



2(i-^V(V)-^V(V)) 

^ _ ^_/^ _ ^^/^ + (2T//i)(/(^+/T) + f{uj^/T)) _ /i 

4 2Ji 



2S (B5a) 
(B5b) 



u.-uj^ + 2T{f{u;+/T) - f{u./T)) = (B5c) 
Here the function / is defined by 

/(x) = -ln[l-e-^] (B6) 
We use eq. ( |B5b| ) to solve for /i. Substituting and rearranging gives 

u;^ + u;_- 2T{f{uj+/T) + /(cu./T)) - i4±i^ = ^ (B7a) 
u;_-co^- 2T{ficojT) - /(CU+/T)) = ^^^^^ (B7b) 

with 

j = 2Ji(l + 8(5,-^)/7r) (B8) 

and 



2nJi{S,-S) 



(B9) 



1 + 8{S, - S)/ti 

These two equations may be easily solved numerically for ujj^ and uj^ by adding the two 
equations to obtain an expression for ujj^ in terms of uj_ and T, and then substituting that 
into one of the two equations to get a single equation for uo-. At low T and sufficiently close 
to the phase boundary the equations have two solutions, one with uj^ < uj_ and one with 
u;_|_ = co'_; above a critical temperature the two solutions merge. By substituting the results 
in to eq. (|A9|) we have verified that where the solution with u;+ < uj- exists it has a lower 



energy than the solution with uo^ = uj_. The results displayed in fig. 3 were obtained in this 



manner. 
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APPENDIX C: SUSCEPTIBILITIES AND RELAXATION RATES 



We begin with the dynamic susceptibhties, which we obtain by computing the hnear 
response of the system to an externally applied magnetic field h\"'\ The Schwinger boson 
formalism is rotationally invariant. We therefore compute only the response to a field in the 
2;-direction. Thus we add to the Hamiltonian a term 

AH = ^hf''^-St^ (CI) 

ia 

After using eq. (|2.2|), eq. (|2.3|), and eq. (^) this becomes 



2 



Ai/ = E '" 01 + : ^ of (C2) 



with 



O^" = Ecosh(0fe+,+P + 0,)(4 

kafS 

+ sinh(4+g+p + 0k){si+g+pa(raf^slkp + h.c. + a ^ s) (C3a) 
'^g = E cosh(^fe+g - 9k){sl^g^(r'^pak+Pi3 + al^^^a'^^Sk+p/s + h.c.) 

kal3 

+ sinh(4+g - ^fc)(s^+^„a^^alfc_p^ + 4+g„<^slfc_^ + h.c.) (C3b) 
The only non-zero correlation functions are 



b{ujk) - b{Uk+q+p) 



4Ecosh2(^fc+g+P + 0^ 



u -Uk + ujk+q+p + ie 



(cufc + Uk+q+p)'^ -{uj + ie)"' 



and 



Xl'iio) = dte"'+"" < [O^{t),O'l,(0)] > 

Y (^^fc + OJk+q)^ -{0J + le)^ 
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We are interested in low energy phenomena; this imphes that k and k + q + P are near or 
P. In this case we may approximate: 

cosh(^fc) = {fi/2ukY/\l + (C5) 

smh{e,) = sgn{^k)if^/2uJkV^^{l - uJkM (C6) 
Then near g = we have 



Xq [^) - 2^ 



k 



^k^k+q LO — U}k + iUk+q ~ 



^ „ (cUfc - Uk+qY (1 + bjUk) + h{uJk+q)){MJk + ^fc+g) 
Y ^k^k+q {uJk + ^k+qY -{UJ + 



aa( X _ (^^fc + ^k+q+pY b{uJk) - b{uJk+q+p) 

Xq [^) — 



^k^k+q+P UJ — LUk + ^k+q+P — "^C 
^ {uJk - UJk+q+pf {I + b{^k) +b{^k+q+p)){i^k + i^k+q+p) .^^N 
^ UkUJk+q+P [UJk + UJk+q+Pj " (^^ + ^(^) 



while near q = P 



k ^k^k+q 



b{uJk) - b{uJk+q) _^ (1 + b{i^k) + b{Uk+q)){u}k + ^k+q) 



UJ -UJk + UJk+q - le [uJk + UJk+q) ^ + 

/i^ [ b{uJk) - b{uJk+q+p) , {l + b{ujk) + b{uJk+q+p)){uJk + i^k^ 



(C9) 



k ^k^k+q+P 



q+P) 



io -ujk + ujk+q+p - «e {uJk + uJk+q+pY - {uj + it)^ 

(CIO) 



In all of these formulae there are important contributions from k near and k near P. 

We now consider NMR. To derive nuclear relaxation rates we take the standard [ |^ 
hyperfine Hamiltonians describing how nuclei are coupled to the electronic spins, write the 
spins in terms of bosons, and then compute the appropriate boson correlation functions. We 
assume throughout that T ^ Ji. 

We begin with the yttrium. An yttrium nucleus sits halfway between two nearest- 
neighbor Cu02 planes and above the center of a plaquette formed by four Cu atoms. We 
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denote the hyperfine coupling to one spin by D. Thus we write the hyperfine Hamiltonian 
for yttrium, 

Hlf^D Y: -^f' (Cll) 

a,i=1..4 

After performing the transformations of section II and Appendix A and retaining only those 
terms capable of giving dissipation at NMR frequencies we have 

Hhf = D J2 9i<l)cosh{9k+g-9k){sl^g^aal3ak+p0 + h.c.) (C12) 

ki,q,a,/3 

with 

\g{q)\ = 4cos{qj2)cos{qy/2) (C13) 

Here we have omitted an unimportant phase factor in g. We now calculate the relaxation 
rate in the usual way, from 

hm - / e^(-+-)* < K/t), i/^/O)] > (C14) 
<^-^o UJ Jo J J 

finding 

-i- = ^ V lg(g)Pcosh^(^^ -^^^.)'^(^^-^^+.) = 2vrD^ hm i V \g{q)\\-" M 

(C15) 

We now consider the planar oxygen. Each oxygen is located in a Cu02 plane and is in the 
center of a bond connecting two Cu sites. Thus 



i=l,2 



Expressing the spins in terms of bosons as was done for yttrium gives: 

C 

^hf = ^ E /(?)( cosh [Ok + dk+q+p){sl^aaj3Sk+q+Pl5 + al+Pa^^ 
k,q 

+ cosh {Ok - Ok+q)s\^(Tapak+q+pp + ^.c.) (C17) 

with 
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|/(A;)|=2cos (I) 



Proceeding as we did with yttrium yields eq. ( |3.10| ) of the text. To evaluate this it is 
convenient to consider q near and q near P separately. The case of q near goes through 
just as for yttrium, except that the square of the form factor is 4, not 16, and one must add 
^aa" _ q near P it is convenient to write g = P + ^2 — fci and to sum over ki and /c2 • The 
integrals have contributions from ki, k2 near and P. The form factor becomes 

\f{h - = {kl - k^,f = {kl + kl)/2 (C18) 

where in the second equality we have done the angular integral and k stands for either k or 
{k — P) as appropriate. Now we have, from eqs. (|A9| ) and ( |2.?| ), 



k^ = A{l-^u)=2{ul-uJ+^_)/^^'' (C19) 

where the gap is uJj^ for k near and uj^ for k near P. Putting this into eq. ( p.lOp yields 
eq. ( |3.11| ). Finally, we consider the Cu relaxation rate. A Cu nuclear moment is believed 
to be coupled to the spin on the same site, via a hyperfine coupling A, and to the spins on 
the four nearest neighbor sites in the same plane, via a hyperfine coupling B. Thus 



i/f; = As','^ + by: sl'^ (C20) 

i=1..4 

After transforming to the boson representation we have 

^hf" = IY.I^- 457(g)] ( cosh {9k + 0k+cj+p){sl^(r^,3Sk+q+P(3 + aLf^a/3«fc+g+P/3) 

+ cosh {Ok - ek+q)sl^g^aaf3ak^g+pp + h.c.) (C21) 



^ kq 



Again we may construct the relaxation rate. It is given by eq. (|3.12| ) and may be simply 
evaluated because the dominant contribution is at q near P. 
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FIGURES 

FIG. 1. Model system considered in this paper: two square arrays of spins with in-plane cou- 
pling Ji and between-plane coupling J2. 

FIG. 2. r = Phase diagram of eq. (^]^) as described in Appendix B from the Schwinger 
boson mean field theory. 

FIG. 3. Temperature dependence (in units of J2) of the Schwinger boson gaps u;+ and 
defined in eqs. ( p. 10 ) and ( 2.11 ) and computed as described in Appendix B, for parameters such 



that at r = the model is at the phase boundary for small J2 and Sc — S. The temperature 
T* ~ J2 separates the low-T coupled planes regime from a high-T decoupled planes regime; in the 
mean field theory there is a second order phase transition at T*; we believe fluctuations would 
convert this to a smooth crossover. 

FIG. 4. Temperature dependence of uniform susceptibility and copper, oxygen and yttrium 
relaxation rates calculated from mean field theory for parameters used in constructing fig. 3. The 
temperature T* at which the model crosses over from the decoupled-planes to the coupled-planes 
critical regimes is indicated. 

FIG. 5. Frequency dependence of antisymmetric and symmetric susceptibilities calculated 
from mean field theory for several different temperatures. 

FIG. 6. Different regimes of behavior of scaling theory and temperature dependence of relax- 
ation rates in each regime. 
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